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A class of three-dimensional models which satisfy supersymmetric intertwin- 
ing relations with the simplest - oscillator-like - variant of shape invariance 
is constructed. It is proved that the models are not amenable to conven- 
tional separation of variables for the complex potentials, but their spectra 
are real and equidistant (like for isotropic harmonic oscillator). The special 
case of such potential with quadratic interaction is solved completely. The 
Hamiltonian of the system is non-diagonalizable, and its wave functions and 
associated functions are built analytically. The symmetry properties of the 
model and degeneracy of energy levels are studied. 
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1. Introduction. 

The idea of shape invariance was proposed [lj in the framework of one- dimensional Su- 
persymmetrical Quantum Mechanics [2], [3] providing a very elegant method of complete 
solution of a class of one-dimensional quantum models. Although almost all these models 
were already solved, this approach gave a new (purely algebraic) look on the class of exactly 



a E-mail: m.bardaveli@gmail.com 
b E-mail: cannata@bo.infn.it 
c E-mail: m.ioffe@pobox.spbu.ru 
d E-mail: cutaisi@yahoo.com 



1 



solvable models. The Hamiltonian H^ 1 ' (x; a) depending on a space variable x and a con- 
stant parameter 7 is said to be shape invariant iff: 1) it satisfies supersymmetric intertwining 
relations: 

H {1 \x; j)q(x) = q(x)H {2 \x; 7), (1) 
with some intertwining operators q, and 2) two Hamiltonians in (JT|) are related as follows: 

H^\x^) = H^\x- 1 ) + R{ 1 ) ) (2) 

where Ri^y) does not depend on coordinate x, and the parameters 7 and 7 are related: 
7 = 7(7). Eq.([2]) just means that intertwined Hamiltonians have the same shape, and they 
differ only a little bit due to a change of parameter. The general form of shape invariance 
was investigated in [I], where all shape invariant potentials with additive dependence of 7 
on 7 were found. 

The method of shape invariance was generalized also to the two-dimensional Schrodinger 
operators not amenable to conventional separation of variables in [6]. It was used in 
formulation of two new - supersymmetric - methods of separation of variables. One of them 
[5], [6] was useful in the study of several quasi-exactly solvable models where a part of 
spectrum and corresponding wave functions was found analytically for arbitrary values of 
parameters of the models. The second one [TJ, [8] provided complete solution of these models, 
but only for specific values of parameters. The general form of two-dimensional systems with 
additive shape invariance was found in [9j. In particular, the new shape invariant potential 
(two-dimensional generalization of Scarf II) was built, and its spectrum was constructed [10] 
by the method of supersymmetric separation of variables mentioned above. 

Of special interest are the models with the simplest kind of shape invariance with 

7 = 7; R(j) = 2X. (3) 

Just this property with intertwining operators q of first order in derivatives is satisfied by 
one-dimensional harmonic oscillator leading to the well known equidistant spectrum with the 
spacing 2A. Much more interesting models with the same properties ([3]) but with higher order 
intertwining operators were studied in [11]. It was shown that depending on the structure 
of zero modes of higher order intertwining operators q and q^ the spectrum of the model 
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consists of finite or infinite sequences of equidistant levels with the spacing 2A. The systems 
with shape invariance of the form ([3]) are naturally called "oscillator-like" shape invariant 
(sometimes also called "self-isospectral"). 

In two-dimensional space, the oscillator-like shape invariance (|3]) with first order differ- 
ential operators q is fulfilled in a nontrivial model with quadratic interaction with such a 
choice of coupling constants that standard separation of variables is impossible (see [T2]). 
The Hamiltonian of this model is non-Hermitian and non-diagonalizable, but its spectrum 
is the infinite equidistant sequence of real energy levels with the same spacing 2A. Moreover, 
the analytical expressions for all wave functions and corresponding associated functions (non- 
diagonalizable Hamiltonian !) were obtained. Quite unexpectedly the same properties of the 
spectrum are present for two-dimensional anharmonic oscillator of some special form [13] as 
well (see also [H]). 

It is known that in Hermitian Supersymmetric Quantum Mechanis with first order su- 
percharges, the generalization [15] to higher dimensionality d > 2 leads [H] to systems with 
separation of variables. But the situation is less restrictive and more interesting for the case 
of complex potentials - the first order operators q may intertwine nontrivial Hamiltonians 
not amenable to separation of variables. Just this opportunity was realized in [12] and [13] 
for two-dimensional models mentioned above. 

In the present paper a class of three-dimensional models with oscillator-like shape invari- 
ance is studied. It is proved that the models are not amenable to conventional separation of 
variables for the complex potentials (Section 2), but their spectra are real and equidistant. 
The special case of such potential with quadratic interaction (Section 3) is solved completely. 
The symmetry properties of the model and degeneracy of energy levels are studied in Sec- 
tion 4. The Hamiltonian of the system is non-diagonalizable, and its wave functions and 
associated functions are built analytically in Sections 5 and 6. 
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2. Three-dimensional oscillator-like shape invariance: general solution. 

The intertwining relations for Hamiltonian with oscillator like shape invariance and with 
supercharges of first order in derivatives, 

HA+ = A + {H + 2A); (4) 
A + = <Xi(x)di + a(x)\ H = —df + V(x); di = d/dxf, i = 1,2,3. 

are equivalent to the system of differential equations for the coefficient functions di(x), a(x) 
and potential V(x) : 

d k di(x) + dia k (x) = 0; (5) 
A (3) ai(f) + 2d i a{x) = -2Aoi(f); (6) 
A (3) a(x) + ai(x)diV{x) = -2Xa{x), (7) 

where A^ 3 ^ is the three-dimensional Laplacian. It follows from ((5]) that 

d1a k (x) = d1a k {x) = d1a k (x) = 0, k = 1, 2, 3, 

and therefore, from (J6j): 

did(x) = — Aaj(x). 

These equations are self-consistent if: 

dia k (x) = d k di(x), i ^ k. 

Thus, the solution of the system (JSJ) - © is: 

di = const; d{x) = —XdiXi, 

where summation over repeated indices is implied, and possible constant in the r.h.s. can 
be eliminated by a suitable shift of coordinates. The remaining equation (J7J) gives: 

didiV(x) = 2X 2 diXi. 

Looking for potential V(x) in the form V(x) = X 2 x 2 + v(x), one has: 

didiv(x) = 0. (8) 



Obviously, at least two of the constants Oj do not vanish, since otherwise variables in Hamil- 
tonian are separated (systems with separation of variables are not the subject of this paper). 
For defmiteness, let us take a\a 2 7^ 0, leading to the general solution of (jHJ): 

v(x) = u(a 2 Xi — aix 2 , a 3 Xi — a\X 3 ) 

and the initial Hamiltonian of the form: 

H = — A*- 3 ) + X 2 x 2 + u(a 2 Xi — dix 2 , a 3 xi — aix 3 ). (9) 

The form (j2J) suggests the opportunity to separate variables in terms of 

yi = a 2 x 1 - a x x 2 , y 2 = a 3 x 1 - a x x 3 , y 3 = hx h (10) 

where b{ are constants. Such transformation of variables is possible if the corresponding 
Jacobian does not vanish: 

D = bia\ + b 2 aia 2 + b 3 ai<x 3 ^ 0. 

Rewriting the Hamiltonian in terms of : 

H(y) = -(a? + a 2 2 )d 2 yi - (a\ + a 2 3 )d 2 2 - b 2 d 2 3 - 2a 2 a 3 d Vl d y2 - 2(a 2 b 1 - aib 2 )d yi d y3 - 
-2(0361 - 0163)^,9^, + \ 2 x 2 k {y) + $(yi, 2/2), (11) 

one obtains that separation is possible in the case of 



b 2 = — } 63 = ^, (12) 



0261 , a 3 6i 

o 2 = 

so that D = 61 a?. Since if the constants a, are real, the Jacobian D does not vanish, the 
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interesting situation not amenable to separation of variables may appear only for a 2 = 0, 
i.e. for complex values of Oj (61 can not be chosen vanishing due to ( Tl2l) . ( fTQl) ), This is a 
reason why only the systems with complex potentials (of = 0) will be considered below (see 
the review papers [17], [IS]). 

One can check that for Oj, a(x) defined above, the intertwining relation (j4j) is fulfilled 
simultaneously with its partner relation, i.e.: 

[H,A ± ] = ±2XA ± , (13) 



where A ± are defined as: 



A ± = ctidi =F XctiXi (14) 

(one has to notice that the operators A^ 1 are not necessarily mutually Hermitian conjugated). 
It is clear from f lT3|) that increase and decrease the energy by 2A, correspondingly If the 
ground state of the model (with zero energy) exists, it is defined by two equations: 

A-^ = {ciidi + XaiX^o = 0, (15) 
H^ = E ^ ; E = 0. (16) 

It is convenient to look for the solution of (j!5p in the form: 

= exp(-Af 2 /2)^ (x). (17) 

After substitution into ( !T5|) . the new function i[) satisfies 

aiditpo^x) = 0, 

and its solution depends on two combinations of coordinates: 

ipo(x) = ipo(a2Xi — aiX2, a 3 xi — CL1X3); \l/o = exp(—\x 2 /2)ip (a2Xi — CL1X2, 03X1 — 01X3). 
According to ( fTBl . ipofaxi — a.1%2, 0,3X1 — 01X3) has to satisfy: 

- A (3 Vo + ^Xidiipo + u*Po = -3A^ , (18) 

which is difficult to solve for ipo with an arbitrary function u. We will act vice versa, expressing 
u in terms of given ipQ. For this purpose, let us write ipo in the exponential form: 

ip = expW(y 1 ,y 2 ), 

and then from ( Fl8l) . we find that: 

u = (diW(x)) 2 + A {3) W(x) - 2Xx i d i W(x) - 3A. (19) 
Thus, the Hamiltonian 

H = -A® + X 2 x 2 + (d.Wix)) 2 + A (3) W{x) - 2Xx l {d l W) - 3A (20) 
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with arbitrary function W has the ground state: 

\I/o = exp(— Ax 2 /2 + W(a2X\ — aix 2 , ^3X1 — 01X3)) (21) 

with zero energy. The excited states with energies E n = 2Xn for the considered case of = 
can be built by the "creation" operator A + : 

* n = (A + ) n ^o = (-2X) n (a iXi ) n %. (22) 
From the definition (fT4j) one has: 

[A + ,A-] = 2\a 2 = 0, (23) 

and therefore, not only the ground state wave function but also all excited wave functions 
are annihilated by A~ : 

A-^ n = 0. (24) 

The full analysis of the general model f!20p seems to be too difficult a task. In particular, 
there is no guarantee that the states ff22l exhaust the entire set of excited states. For this 
reason, in the next Section a particular case of the model will be studied in more details. 

3. The particular case. 

Let us consider here the particular case of parameters of the model: 03 = 0, a\ — 1, 02 = — i, 
and the simple form of function W which leads to the quadratic interaction: 

W(a 2 x 1 - atx 2 , x 3 ) = g(x x - ix 2 )x 3 = gzx 3 , 

where z = x\ + ix 2 z = X\ — ix%. Then, from (fTI]) . fT20l) . f[2Tj) one has: 

H=-8$ + \ 2 x\ + g 2 z 2 - 4:\gzx 3 - 3A, (25) 
A ± = 2d z =f \z, (26) 
^0 = exp(-Af 2 /2 + gzx 3 ), (27) 

where the mutual independence of z, z provides (1231) . and for A > \g\ the wave function is 
exponentially decreasing at infinity. 
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Actually, this system reproduces exactly one of the three-dimensional models of a recent 
paper [19], if in Eq.(35) of that paper b = 2A, b 3 = 2g, x 2 — > — x 2 and Q ± -H- Q T . Thus the 
same Hamiltonian (1251) is intertwined by supercharges Q ± of second order in derivatives: 

[H,Q ± ] = ±AXQ ± , (28) 

where 

q- = d\ + d\ - dl + ca + b, 

g+ = dl + <9 2 2 - dl - ca + b- d^, 

C 1 = 2(\x 1 - gx 3 ), 

C 2 = -2(Ax 2 + igx 3 ), 

C 3 = 2(-\x 3 + gz), 

B = \ 2 (zz-x 2 3 )-g 2 z 2 + X, 

CAi = 2X(zd z + zd- z - x 3 d 3 ) - Agx 3 d z + 2gzd 3 . 

The intertwining relations (I28I) also have the form of oscillator-like shape invariance of 
Hamiltonian H, but with the double spacing between energy levels. The minimal energy state 
corresponds to the zero mode of Q~, and it can be easily found. Indeed, direct calculation 
shows that: 

Q- = * Q {dl + d 2 2 -di)%\ 
q+ = q--2\- 2CA, 

therefore, the same ground state \l/o given in (1271) is annihilated not only by A~ , but also 
by Q~. One can check that other possible zero modes of Q~ can not be simultaneously the 
wave functions of H. 

It is clear that both (A + ) 2 ty and Q + ^o are the wave functions of the Hamiltonian H 
with the same energy 4A. The calculation shows that 

Q + ^o = -2(A^ + Cid&o) = -2A[1 + 2A(x 2 - zz)]V - 8\g 2 z 2 V , 

i.e. there is degeneracy at the second excited level 4A with two independent wave functions: 

^ 2 = (A + ) 2 ^ ~ z 2 ^ 0} $ 2 = [1 + 2A(x 2 - zz)]V . 
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Analogously, higher wave functions of the form Q + ty n (x) can be calculated, using expression 
( 122|) . which, in our case of a\ = 1, a 2 = —i, a 3 = 0, means \l/ n ~ (z) n ^ : 

Q + z n ^ = (Q+ -2X- 2CA)z n ^ = -2Xz n V - 2z n C i d l ^ - 2^ CidiZ n = 
= -2X[2n + 1 + 2X(x 2 3 - zz)]z n ^ ~ 8Xg 2 z 2+n ^ . 

Therefore, for arbitrary n, the function 

^ n = [2n + 1 + 2X(x 2 3 - zz)]z n ^o 

is the wave function of H with energy value 2Xn + 4A, being degenerated with \l/ n+ 2- This is 
the general situation: action of operators Q + leads to degeneracy of levels. Degenerate wave 
functions are: 

^Ux) = (Q + ) k (A + ) n V ; ^ 0n ee E kn = 2X{n + 2k), (29) 

with the multiplicity of level N = k + 1 + [Vi/2]. It is necessary to mention that the order 
of " creation" operators in (129]) is irrelevant due to following commutation relations between 
operators Q ± and : 

[A + , Q~] = 4AA- [A~,Q+] = -AXA + ; (30) 
[A + ,Q + ] = [A-,Q-] = 0. (31) 

One may expect that combinations A ± and Q ± of third order in derivatives will give 
additional degeneracy of some levels, since they also satisfy oscillator-like shape invariant 
intertwining relations obtained from (fl3|) and (128]) : 

[H,A + Q-] = -2XA+Q-, [H,Q-A+] = -2XQ-A+; (32) 
[H, A~Q + ] = 2XA~Q + , [H, Q + A~] = 2XQ + A~; (33) 
[H,A + Q + ] = 6XA+Q+, [H,A~~Q~] = -6XA~Q~. (34) 
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However, additional degeneracy is not generated, since: 

A + Q-* n = A + y (d! + d\ - dl)z n = 0, 

Q-A+t> n = (A + Q- - AXA-)^ n = 0, 

A-Q + ^ n = A-(Q- - 2A - 2CA)^n = 

= -2A-GdiV n = -2(CAA- + 2\A + )V n = -AXA + ^ n , 

Q + A~V n = 0, 

(35) 

i.e. action of these third order operators does not give any new states. 
4. Symmetries. 

The degeneracy of levels which was described in the previous Section indicates the existence 
of some symmetry of Hamiltonian. Indeed, the commutation relations (fT3|) . (|28|) lead imme- 
diately to the symmetry operators commuting with the Hamiltonian (notations will be clear 
below) : 

Ro = A + A- = A-A + ; R 1 = ^[Q + ,Q-]; R 2 = Q + Q-. (36) 

These operators are bilinear in "creation" and "annihilation" operators, and they definitely 
do not change the energy of the state. The first of these operators is of second order in 
derivatives, and the others two are of fourth order, but the order of R\ can be reduced by 
two units, since: 

R x = —AXH + lQgd z d z + 8g\z(gz - \x 3 ) - 12A 2 . 

Thus, a more convenient symmetry operator instead of R x is 

R 1 = 2d z d 3 + \z{gz-\x 3 ). (37) 

Also, the fourth order symmetry R 2 can be replaced by third order symmetry operator which 
is proportional to the commutator [Rq, R 2 ). Relations fl30l . fl3Tl) help to express this operator 
as: 

R 2 = A\zd z Q--(\ + C i d i )(A-) 2 . (38) 
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One can calculate the mutual commutators between the symmetry operators Ro,Ri,R 2 . 
They are: 



[R ,R 1 ] = 0; [R 2 , Ro\ = AXR 2 ; [R 2 , R{\ = -2gR 2 . (39) 

More complicated symmetry operators consist of the products of A ± and Q ± such that 
the energy of state is not changed: 

R 3 = Q+(A-) 2 ; R 4 = Q~(A + ) 2 . (40) 

Using the explicit expressions for A ± , Q^ 1 and commutation relations ( 1301) . ( 13T1) . one can find 
that: 

i?3 — -R4 = 2R 2 "I" 8XRq, 

so that only one of them - for example, R3 - can be taken as independent. Its nonvan- 
ishing commutators with i? , -Ri, R2 do not give any new independent symmetry opera- 
tor. The straightforward calculations prove that higher symmetry operators of the form 
(A+) 2n {Q-) n , (A~) 2n (Q + ) n as well as the products with different order of A and Q also do 
not provide new independent symmetries. 

In the context of conventional - Hermitian - Quantum Mechanics, such situation would 
be called as maximally superintegrable [20] . Indeed, for the three-dimensional system one 
has two mutually commuting symmetry operators Rq, R± (complete integrability) and, in 
addition, two symmetry operators R 2 , R3, which do not commute with them. Commutators 
of this last operator with polynomials of Rq, R\ are not functionally independent from Rq 
and do not provide additional symmetries. 

Now, it is instructive to show the connection of symmetries with degeneracy of energy 
levels. Using commutation relations between operators A ± , Q ± and definitions of R$, R±, R 2 , 
one finds the action of the symmetry operators onto wave functions ^kn, ( 1291) : 



Ro^kn = -4AA;* (fc _i) (n+2) ; (41) 

R^ kn = 2gk^ {k _ 1)(n+2 y, (42) 

R2^kn = 8g 2 k(k-l)^ {k ^ 2Kn+A) -AX 2 k(2n + 3)^ {k . 1){n+2) ; (43) 

R 3 V kn = lQ\ 2 k{k - l)m {k _ 1){n+2) . (44) 

11 



It is evident that all wave functions in ( 14"T|) - ( I44j) correspond to the same level with energy 
Ek n = 2X(n + 2k). The repeated actions of symmetry operators span all wave functions which 
belong to this degenerate level. 



5. Norms, scalar products and non-diagonalizability of Hamiltonian. 

The Hamiltonian (|25|) is not Hermitian, but it obeys the property of pseudo-Hermiticity [21 J : 

•qHrf 1 = # f (45) 

with r\ an invertible Hermitian operator. In the present case, an operator r\ can be chosen as 
r] = P2, the inversion x 2 — > —^2- The theory of one-dimensional pseudo- Hermitian systems 
was developed during last decade [2T] , the pseudo-Hermitian models in two dimensions were 
studied in [23], [12], [13], where the basic formulas can be found. In particular, Quantum 
Mechanics for such systems with unbroken i]T = P{T— symmetry P{T^ n {x) = ^ n (x) is built 
with the new scalar product (7]— product), which in the present case has the form: 

< $ I 77 I $ >=<< * I $ >>= J \J/$, (46) 

instead of the standard J 

With the scalar product fj46|) . the norm of ground state \I/o is: 



J ^ld 3 x = J exp(— \%1 + 2gzx 3 )d 3 x = - J exp(—\zz)dzdz J exp(— \xl + 2gzx^)dx^, 



= \\j\j exp(-A*s + 9 -^-)dzdz = a/RaF- (47) 
The norms of excited states from fl22l) can be also calculated: 



« ^ n I ^ n » = « (A + ) n V I (A + )"^ »= (-1)" « ^0 I #0^0 »= 

= (-l) n « ^0 I ^0 » <W (48) 

In turn, the norms of excited states — {Q + ) k ^n given by f l2"9~j) can be found as well, 
although with more involved calculations using Ri^ n = 0; Q~^ n = and the commutation 
relations between Q ± : 

[Q + ,Q-] = -8\(H + 3\) + 16gR 1 . (49) 
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One can check that required scalar products are: 

where the coefficients in the last step are not important. Analogously, one can check that 
for different values of indices k, n, the scalar products 

<< ^ kn | ^k'ri >>~ $kk'$nri$nO- 

(51) 

Thus, all wave functions ^>kn of the model, excluding the ground state \l/o and the states 
^ko-, have zero norms and are orthogonal to each other. The self-orthogonality of wave 
functions signals that we deal with a non-diagonalizable Hamiltonians. This situation was 
studied recently in one-dimensional [22] and two-dimensional [12], [13] context, where nec- 
essary technical details can be found. Non-diagonalizability means that each wave function 
^kn, n > with zero norm must be accompanied with a set of associated functions which 
participate in the resolution of identity. 

6. Associated functions. 

Referring back to the papers [22] , [12] , [13] for details on the structure of non-diagonalizable 
non-Hermitian Hamiltonians, biorthogonal basis, scalar products between wave functions 
and associated functions and nontrivial form of resolution of unity, let us remind the most 
important relations which must be fulfilled. 

For such systems, each self-orthogonal wave function ^> kn {^) with zero norm (below they 
be denoted as ^kn,o(x)) must be accompanied with a set of p kn — 1 associated functions 
^kn,m, m = 1,2, ...,pkn — 1? where p kn is called the dimension of Jordan cell. This situation 
should be distinguished from the conventional degeneracy of the energy level (here this 
degeneracy also exists). A new notation with an additional index of wave function ^> kn = 
^kn,o is useful, since by definition these functions obey: 

(H - E kn )^kn,rn = ^kn,m-l, TH = 1, 2, p kn - I] (H - -Ejfc n )*fcn,0 = 0, (52) 

where all functions are supposed to be normalizable according to the scalar product ( 1461) . 
Thus, each self-orthogonal eigenfunction ^kn,o is supposed to be accompanied by a set of 
associated functions ^kn,m, m = 1,2, ...,Pk n ~ 1- 
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With these notations, according to the general formalism which was given in detail for 
some one-dimensional models [22] and two-dimensional models [12] . [13], the scalar products 
must be: 

((^kn,m\^k'n',m')) = J ^ kn,m{x)^ k'n' ,m> {x)d 3 X = 5 kk >5 nn/ 5 m ( Pfen _ m /_i) ; Ul' = 0, 1,2, ...,p kn - 1(53) 

The Hamiltonian H is clearly non-diagonal, but block-diagonal. Each block - Jordan cell 
- has dimensionality p kn , which will be shown (see Eq. (l59l) ) to equal: 

Pkn = n + 1. (54) 



The self-orthogonality of wave functions - ^kn,o, n^O - was already demonstrated in ( 1501) . 
( 148]) . To complete the construction of the Jordan cell, it is necessary to find the corresponding 
associated functions with properties (1521 . (1531) above. The procedure is rather similar to that 

in p, m- 

We will start from ^ kn fl with k = 0. Due to relations 

(A-) k (H - E n ) = (H - E n _ k )(A-) k , 

the associated functions \I/on,m for corresponding wave functions 

Xx 2 

*on,o = (A + ) n ^ 00 ,o = c n , ^ exp ( — ^ + ^X 3 ) 
in (l29l) (c n) o are constants) satisfy, as in [12], [13], the equations: 

where a„ im constants. The solution is searched in the form: 

Xx 2 

^On,m = exp( + gZX 3 )tPn,m, 

where the function ip n ^ m has to satisfy: 

P) m , n - n r ^n-m 

u z rn,m u "n,m u n— mfl^ 

Its general solution is: 



TO— 1 

mi ^— ' 



i=0 
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where functions gn,m{z,Xz) can be found from equations: 

(H - E n )^ 0nim = ^on, m -i (55) 

by direct recurrence procedure. An alternative way is to calculate the associated functions 
for the lowest values of n, and to guess their form for higher values. Indeed, one can prove 
by induction that 

Ar 2 2f7 

*on, m ~ (A+)"- m exp(-^ + gzx 3 ){z - -fx 3 ) m ~ (A + ) n ~ m d^ 0fi , m = 0, 1, .., n 

satisfies fl52l) . This expression leads (by straightforward calculations) to the necessary scalar 
products: 

In turn, the associated functions \l/fcn,m for k ^ satisfy: 

(H - E kn )^ kn>m = * fcnjW _i, £ fen = 2X(n + 2k). (56) 
Taking into account that 

(Q + )\H - E 0n ) = (H- E kn )(Q + ) k , (57) 

one can try to identify: 
It is easy to check that 

= « ^On,m | (Q'f^knfi »~« ^On,m | ^0n,0 »~ <W (58) 

i. e. ^kn,o is orthogonal to all ^kn,m with m ^ n. In order to satisfy the necessary con- 
ditions for all scalar products, one has to use the freedom to add the solution ^kn,o of the 
homogeneous equation ( |52l) with arbitrary constant coefficient. This is enough to fulfill the 
conditions: 

*^*^ ^ kn,m | ^A;n,m' ^m(n— m')j 
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which mean that we deal with Jordan cells of matrix dimensionalities pk n — n + 1. Analogous 
calculations show that: 

<< ^kn,m | ^fc'n'.m' >>= 5 nn '5kk'5m(n-m'), (59) 

completing the proof that the Hamiltonian is block-diagonal with Jordan cells of dimension- 
ality p kn = n+1. 

7. Conclusions. 

In the present paper the method of supersymmetric shape invariance was applied to in- 
vestigate a new model. This model has several interesting properties which distinguish it 
from others recently studied [12], [13]. First of all, it is three-dimensional, and although 
the interaction is quadratic, the special choice of coupling constants with one of them pure 
imaginary does not allow to separate variables. Similarly to two-dimensional models of [32] , 
|13j . the Hamiltonian is non-diagonalizable, but now it incorporates nontrivial symmetries. 
Two of the symmetry operators are mutually commuting, while two others are not. Besides 
the block-diagonal form of Hamiltonian with Jordan cells corresponding to wave functions, 
each energy level is degenerated due to symmetry operators which mix the wave functions 
with the same energy. The number of symmetry operators just implies that our model be- 
longs to a class of maximally superintegrable models in conventional Hermitian Quantum 
Mechanics, but we do not know of general statements on possible numbers of symmetries in 
non-Hermitian case. The analysis above may lead, in particular, to some hints concerning 
such general theorems. 
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